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ABSTRACT 

The  effects  of  varying  the  configuration  of  the  Large  Aperture  Seismic 
Array  (LASA)  on  the  detection  signal-to-noise  ratio  observed  in  the  SeismLc 
Data  Analysis  Center  Detection  Processor  are  discussed.  It  is  shown  that  a 
configuration  for  the  LASA  of  thirteen  subarrays  (A,  B,  C,  and  D rings)  with 
16  sensors  per  subarray  produces  an  average  loss  of  0.2+0. 3 db  (one  standard 
Jtviation  of  the  mean)  from  the  current  configuration  of  seventeen  subarrays 
(A,  B,  C,  D,  and  E rings)  using  sixteen  sensors  per  subarray.  Reducing  the 
number  of  subarrays  to  nine  (A,  B,  and  C rings)  with  only  7 elements  per 
subarray  produces  an  average  loss  of  3. 5+0. 3 db.  These  results  are  in 
general  agreement  with  known  results  on  variation  of  noise  reduction  and 
signal  loss  as  a function  of  subarray  and  subarray  sensor  selection. 
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INTRODUCTION 


A series  of  off-line  Detection  Processor  (DP)  experiments  were  performed 
on  eata  from  the  Large  Aperture  Seismic  Array  (LASA)  in  Montana  to  determine 
the  effects  that  the  number  of  subarrays  and  number  of  sensors  per  subarray 
have  on  the  output  of  the  cc — uterized  seismic  analysis  system  at  the  Seismic 
Data  Analysis  Center  (SDAC)  in  Alexandria,  Virginia.  The  outputs  of  the 
system  are  the  seismic  events  listed  in  a Daily  Summary,  the  DP  detections, 
etc.  (See  Dean  et  al.,  1971,  for  a complete  discussion  of  the  SDAC/LASA 
system) . 

The  functional  flow  of  data  from  LASA  through  the  SDAC/LASA  system  is: 

1)  data  acquisition;  2)  detection  processing;  3)  event  processing;  4)  experi- 
mental operations  console  (EOC)  editing;  and  5)  publishing  of  the  Daily 
Summary.  The  experiments  discussed  in  this  report  concern  themselves  only 
with  function  2),  detection  processing,  and  the  resultant  signal-to-noise 
ratios  of  the  detected  signals.  The  detection  parameters  used  are  those 
discussed  by  Chang  (1974)  using  LASA  Beam  Set  1_>3  in  Partition  I and  LASA  Beam 
Set  140  in  Partition  II.  (Experiment  7,  discussed  below  was  an  exception, 
using  LiiS  133  also  in  Partition  II.) 

Earlier  studies  of  the  effects  of  reducing  the  number  of  elements  at 
LASA  (Hartenberger , 1967;  Hartenberger  and  Van  Nostrand,  1970)  have  shown  that 

the  signal-to-noise  ratio  loss  is  less  than  2db  compared  to  that  for  the 

original  525  sensor  array  when  the  number  of  elements  is  reduced  to  119  or  51 

with  minimum  sensor  spacing  of  3km  or  6km  respectively.  All  of  the  data  used 

in  the  earlier  studies  were  prefiltered  (0. 4-3.0  Hz),  were  beamed  to  the  known 
epicentral  locations,  and  were  corrected  for  travel-time  anomalies  (Chiburis, 
1968).  Also,  the  event  set  contained  only  earthquakes  well  above  magnitude  4.7. 

Such  differences  (between  the  present  study  and  these  earlier  ones)  as 
filter  pass-band  (0.9-1. 4 Hz  here)  travel-time  residuals,  and  subarray  and 
array  beam  deployment  may  adversely  affect  an  accurate  comparison  between 
them.  The  comparisons  of  relative  S/N  improvement  between  experiments  in 
this  paper  should,  however,  be  valid,  since  these  parameters  are  held  constant 
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between  experiments.  Care  must  be  taken,  however,  that  non-detection  of  weak 
events  by  small  subsets  of  the  full  LASA  does  not  lead  to  biased  estimates  of 
the  (S/N)  loss.  The  experiments  performed  in  this  report  simulate  the  on-line 
processing  of  a continuous  data  stream  and  include  events  near  the  detection 
threshold  of  the  array  (m^  = 3. 7-4.0) 

For  reference,  the  configuration  of  LASA  is  shown  in  Figure  1,  and  a 
subarray  configuration  in  Figure  2. 

In  order  to  vary  any  of  the  several  DP  parameters  j.n  the  system  for  com- 
parative analyses,  an  off-line  DP  set  of  programs,  written  by  IBM  under 
contracts  F19628-67-C-0198  and  F19628-68-C-0400,  was  usee  to  simulate  the 
on-line  programs.  The  off-line  DP  results  were  then  compared,  detection  by 
detection,  and  the  effects  of  changing  the  array  configuration  were  evaluated 
statistically. 
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Figure  2.  LASA  svbarray  configuration. 


PROCEDURE 


rhe  data  selected  for  analysis  are  those  48  events  appearing  on  the  LASA 
bulletin  for  the  time  period  May  22,  1930Z  to  1972  May  23,  2040Z.  A list  of 
these  events,  taken  from  the  LASA  bulletin  is  given  in  the  Appendix.  This 
particular  time  period  was  selected  because  it  had  been  analyzed  earlier 
(Ahner,  1973)  in  the  EP  at  a lower  signal-to-r oise  threshold  (10  db)  than  is 
normally  used,  and  we  felt  that  there  would  be  no  surprises  in  the  data. 

Of  the  48  events  in  the  Appendix,  7 are  marked  by  an  X and  were  not 
considered  because  their  (S/N)  ratio  was  less  than  14db  on  the  full  array 
bectn  detections.  This  reduction  helped  r inimize  errors  due  to  noise  contami- 
nation of  the  signal.  As  we  shall  see,  however,  far  fewer  events  than  the 
remaining  41  were  involved  in  the  final  calculations.  The  actual  numbers 
range  between  9 and  23.  Many  events  were  missed  because  of  tape  reading  errors. 
These  were  by  no  means  consistant  from  experiment  to  experiment;  and  as  a 
result  in  several  runs  which  should  have  been  identical,  different  numbers  of 
evenLC  were  detected.  Whenever  a tape  reading  error  was  encountered,  the 
computation  of  the  long-term  noise  average  (LTA)  was  disturbed.  This  also 
would  introduce  variations  of  detections  from  run  to  run.  To  minimize  this 
effect;  no  events  which  occurred  within  5 minutes  of  a tape  reading  error  were 
considered.  The  LTA  computations  also  had  to  be  restarted  each  time  one  of 
tha  more  than  30  tapes  covering  the  time  period  came  on  line.  In  the  runs 
with  only  a few  sensors  operational  several  events  would  of  course  simply  be 
missed  because  of  the  lowered  threshold. 

Another  reason  for  differing  results  between  runs  is  that  the  LTA  is 
updated  only  every  third  0.6  second  tinv-step.  For  some  start  times  (after 
a faulty  tape  read  for  example)  the  LTA  may  be  contaminated  by  the  signal. 
Contamination  of  the  LTA  is  in  evidence  in  Figure  3 which  shows  the  STA/LTA 
as  a function  of  the  STA  for  the  37  events  detected  by  the  full  array  in 
off-line  Experiment  1 (to  be  discussed  below).  The  event  marked  by  X shows 
a severe  contamination.  ' .oughout  this  report  events  were  removed  from  the 
plots  and  from  calculation  of  mean  S/N  differences  if  the  S/N  difference  with 
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respect  to  a standard  run  deviated  from  the  mean  by  more  than  two  standard 
deviations  of  an  indi  iiual  observation.  This  criterion  lead  to  the  loss 
of  2-4  events  per  experiment.  No  bias  should  result  from  this  procedure 
since  the  error  is  equally  likely  to  occur  for  either  of  two  different 
experiments  being  compared. 

It  is  noticable  that  the  slote  of  the  data  points  in  Figure  3 is  less 
than  1.0.  This  presumably  arises  because  the  LTA  will  be  contaminated  more 
if  the  signal  under  consideration  is  large. 

The  series  of  experiments  selected  to  demonstrate  the  effect  of  subarray 
and  sensor  reduction  is  listed  in  Table  I.  In  normal  operations  the  on-line 
DP  uses  two  different  array  configurations.  The.  first,  called  Partition  I, 
uses  seventeen  subarrays  from  the  A,  3,  C,  D,  and  E rings,  and  the  second. 
Partition  II,  uses  nine  subarrays  from  the  A,  B,  and  C rings.  Both 
partitions  utilize  sixteen  sensors  per  subarray. 
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PREI.IMINARY  ANALYSIS 


An  initial  observation  of  the  effects  cf  array  size  can  be  obtained, 
without  reporting  to  the  experiments  listed  in  Table  I,  by  using  only  the 
results  of  the  or-line  DP  and  comparing  Partition  I with  Partition  II. 

Figure  4a  shows  the  on-line  signal-to-noise  ratio  from  Partition  II  (ordinate) 
compared  to  Partition  I (abscissa).  As  can  be  seen,  the  (A,  B,  C)  144-element 
array  of  Partition  II  has,  on  the  average,  a signal-to-noise  ratio  2. 0+0. 2 
(one  standard  deviation  of  the  mean,  db  less  than  the  (A-E)  272-element 
Partition  I. 

The  array  configuration  in  Partition  I for  Experiment  1 is  identical  to 
that  of  Partition  I of  the  on-line  array.  In  Figure  4b  we  compare  off-line 
to  on-line  results.  We  expect  to  see  little  difference  if  the  off-line  DP 
package  is  performing  correctly.  We  observe  no  difference,  and  in  fact  the 
average  loss  in  Experiment  1 was  2.0db.  Hence,  it  is  assumed  that  the  off-line 
DP  programs  are  adequate  for  simulating  the  on-line  DP  programs.  All  subsequent 
reduced  array  experiments  are  compared  to  the  Partition  I off-line  simulation 
in  Experiment  1. 
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SIGNAL  TO-NOISE  RATIO  (17.16) 

Figure  4a.  Comparison  of  the  two  detection  beamsets  in  the  on-line 
system  (Partition  II  vs  Partition  I).  Intercept  of  the  line  with 

unit  slope  is  determined  by  the  average  difference  in  log  (S/N)  between 
partitions. 
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RESULTS  OF  EXPERIMENTS 


The  a*. ray  configuration  for  Experiment  2 retains  seventeen  subarrays 
tor  Partition  I and  nine  subarrays  for  Partition  II,  but  the  number  of  sensors 
per  subarrjy  is  reduced  from  sixteen  to  ten  by  dropping  the  inner  two  rings 
of  sensors  from  each  subarray  (see  Figure  2).  The  results  are  shown  in 
Figure  5.  The  signal-to-noise  ratio  loss  relative  to  Experiment  1,  due  to 
the  sensor  reduction,  averages  O.gdt  for  Partition  I and  2.3db  for  Partition  II. 


The  configurations  for  Experiment  3 are  seventeen  and  nine  subarrays  with 
seven  sensors  per  subarray,  obtained  by  dropping  the  inner  three  rings  in 
each  subarray.  The  results  are  shown  in  Figure  6.  The  relative  loss  for 
Partition  I averages  i.8db,  and  for  Partition  II,  3.5db. 


For  Experiment  4,  the  array  aperture  is  reduced  to  one  half  by  dropping 
the  E-ring  subarrays  (leaving  a total  of  thirteen  subarrays),  but  the  number 
of  sensors  per  subarray  is  returned  to  sixteen.  The  results  are  shown  in 
Figure  7 where  the  Partition  I loss  averages  only  0.2db  and  the  Partition  II 
loss  averages  2.6db.  The  significant  result  of  Experiment  4 is  that  the 
signal-to-noise  ratio  losses  are  so  small  when  the  E-ring  is  dropped  from 
the  array.  The  reason  for  this  result  is  that  any  loss  due  to  a higher  noise 
background  with  the  smaller  number  of  sensors  is  off-set  by  better  signal 
correlation  between  subarrays  which  are  closer  together  (Hartenberger  and 
Van  Nostrand,  1972). 


The  configuration  for  Experiment  5 is  similar  to  Experiment  4,  except 
that  the  number  of  sensors  per  subarray  is  reduced  from  sixteen  to  thirteen 
(inner  ring  dropped  from  each  subarray).  The  results,  shown  in  Figure  8, 
indicate  a Partition  I average  loss  of  0.4db  and  a Partition  II  average  loss 


of  2.5db. 


Reducing  the  number  of  sensors  per  subarray  to  ten  and  using  thirteen 
subarrays  (Experiment  6),  yields  the  results  shown  in  Figure  9.  The 
Partition  I loss  averages  0.8db  and  the  Partition  II  loss  averages  3.1db. 
Experiment  7,  with  the  lowest  total  number  of  sensors  in  the  entire  set  of 
experiments,  has  a configuration  of  thirteen  subarrays  (Partition  I)  and 
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re  8.  Experiment  5 result 


nine  subarrays  (Partition  II),  both  with  seven  sensors  per  subarray.  However, 
the  nine  subarrays  in  Partition  II  for  Experiment  7 are  composed  of  the  A,  C, 
and  D subarrays,  instead  of  the  usual  A,  B,  and  C subarrays.  The  reason  for 
this  change  is  to  note  the  contributions  to  the  si gnal-to-noise  ratio  made  by 
the  B-ring  subarrays  which  are  close  to  the  A-  and  C-ring  subarrays  (see 
Figure  1).  Since  LBS  140  was  designed  only  for  the  low  aperture  ABC  ring 
LASA,  LBS  133  was  inserted  in  Partition  II  for  this  Experiment.  The  results, 
shown  in  Figure  10,  indicate  a Partition  I average  loss  of  i.8db  and  a 
Partition  II  average  loss  of  ll.ldb.  This  latter  figure  is  far  out  of  line 
with  what  is  expected  by  comparison  with  Partition  II  of  Experiment  3 which 
had  a loss  of  3.5db  for  the  A,  B,  and  C rings.  We  initially  thought  that 
the  trouble  might  be  that  the  direct  sum  of  the  individual  traces  was  so  small 
that  truncation  error  in  the  computer  was  the  cause  of  the  problem.  Therefore 
we  repeated  the  run  with  all  of  the  raw  data  scaled  up  by  a factor  of  4. 

While  15  events  were  detected  in  this  case,  as  compared  to  9 in  the  previous 
one;  the  events  in  common  had  an  average  mean  log  S/N  value  within  .05  units 
of  each  other.  Thus  there  appeared  to  be  no  errors  of  computation.  Close 
examination  of  the  detection  beam  numbers  for  Partition  II  showed  that  approxi- 
mately 80%  of  the  detections  were  not  on  the  same  beam  number  as  the  detections 
in  Partition  I of  Experiment  1;  whereas  for  Experiment  3,  80%  of  the  detections 
were  on  the  same  beam  number.  This  fact  suggests  a large  number  of  side- lobe 
detections.  This  would  be  consistant  with  the  fact  that  elimination  of  the 
B ring  subarrays  will  sharpen  the  main  lobe  of  the  full  array  at  the  expense 
of  greatly  increasing  the  side  lobe  amplitudes.  Detections  on  the  side  lobes 
will,  of  course,  have  low  signal-to-noise  ratios.  Averaging  over  all  detected 
events  in  both  repetitions  of  Experiment  7 we  obtain  a loss  of  6.2db.  However, 
because  of  the  side  lobe  problem  discussed  above,  these  results  should  not  be 
compared  to  theoretical  predictions.  This  problem  should,  however,  be  kept 
in  mind  in  the  design  of  new  arrays;  the  lesson  is  that  a compact  array  design 
will  make  the  task  of  automatic  detection  simpler. 

Examination  of  Figures  5-11  reveals  the  effect  of  signal  contaminated 
noise.  At  high  S/N  levels  the  points  tend  to  lie  above  the  line.  Since  the 
noise  level  is  lower  for  the  standard  (Partition  I,  Experiment  1),  it 
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exper’.nces  a greater  percentage  increase  due  to  signal  contamination;  and 
this  results  in  a greater  reduction  of  S/N.  Quantitatively  the  effect  is 
small;  omitting  these  points  would  increase  the  loss  by  ~0.3db  in  the  most 
severe  case. 

The  results  of  the  seven  experiments  are  summarized  in  Table  II.  We 
also  give  the  predicted  noise  loss  in  the  pass-band  0. 8-2.0  Hz  by  using  the 
estimation  techniques  of  Blaruford  and  Clark  (1971)  at  the  subarray  level, 
and  assuming  /n  noise  reduction  between  subarrays.  The  difference  between 
these  predictions  and  the  observed  loss  may  be  called  the  apparent  signal 
loss.  (Of  course  we  expect  some  error  here  since  the  actual  noise  and  signal 
pass-band  is  0.9-1. 4 Hz.  Examination  of  Figure  6 in  Blandford  and  Clark, 
together  with  the  realization  that  the  noise  loss  is  controlled  by  the  low- 
frequency  end  of  the  spectrum  and  that  only  relative  loss  is  of  concern  in 
this  analysis  leads  to  the  conclusion  that  the  error  is  on  the  order  of  0.2db.) 

A rough  confidence  interval  on  the  observed  loss  values  may  be  obtained 
by  noting  that  the  average  standard  deviation  of  a signal  loss  estimate  from 
the  mean  of  an  experiment,  averaged  over  Experiments  1-6,  is  1.2db.  Since 
the  average  number  of  observations  per  experiment  is  15,  the  average  standard 
deviation  of  the  mean  would  be  about  0.3db;  and  a 95%  confidence  interval 
would  be  about  0.6db. 

Comparing  Row  1 to  Row  4 in  Table  II  we  see  that  the  apparent  signal  loss 
from  dropping  the  E ring  is  -l.Odb.  That  is,  there  is  less  signal  loss:  as 
would  be  expected.  This  result  is  in  agreement  with  the  reduction  of  signal 
loss  of  0.9db  for  a Ryuku  event  in  the  pass-band  0. 7-2.0  Hz  found  for  dropping 
the  E ring  by  Hartenberger  and  Van  Nostrand  (1972). 

As  the  number  of  sensors  per  subarray  decreases  however  (rows  2 and  3 
compared  to  1;  rows  5,  6,  and  7 compared  to  4)  the  signal  loss  becomes  greater. 
This  is  due,  we  believe,  to  the  fact  that  the  travel-time  residuals  are  set 
for  the  center  sensor  of  the  array.  As  the  number  of  sensors  per  array  is 
reduced,  the  average  contribution  of  sensors  near  the  center  element  is 
reduced;  thus  increasing  the  signal  loss.  For  9 subarrays  this  pattern  is 
not  so  apparent*  we  see  only  that  the  loss  in  row  13  for  7 sensors  per  subarray 
is  less  negative  than  it  is  for  10,  13,  or  16  elements  per  subarray. 
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The  above  discussion  suggests  that  the  observed  loss  figures  have  their 
roots  in  geophysically  real  effects,  and  that  they  can,  therefore,  b_  trusted 
as  the  best  available  estimates  of  the  signal-to-noise  ratio  loss  we  may 
expect  from  such  array  modifications. 

In  Figure  12  we  have  plotted  the  observed  loss  as  a function  of  number 
of  sensors  per  subarray  for  different  numbers  of  si'barrays.  We  see  that 
dropping  the  E ring  results  in  a S/N  loss  of  less  thar  0.2db;  while  dropping 
the  D and  E ring  and  reducing  the  number  of  sensors  from  16  to  7 results  in 
a loss  of  3.5db.  We  have  seen  that  the  standard  deviation  of  these  estimates 
is  approximately  0.3db. 
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CONCLUSIONS 

The  results  obtained  from  the  present  study,  in  which  several  experiments 
were  performed  on  LASA  data  with  reduced  array  configurations,  lead  to  the 
following  conclusions: 

1*  The  detection  signal— to— noise  ratio  for  a 16  element  subarray  D— ring 
size  LASA  is  within  0.2+0.3db  of  that  for  an  E-ring  sized  LASA. 

2.  Further  reducing  the  array  to  an  aperture  of  the  C-ring  produces 
losses  of  2.3j^).3db. 

3.  Eliminating,  within  a subarray,  the  six  sensors  nearest  the  center 
produces  less  than  one  db  of  loss  in  the  detection  signal— to— noise  ratio. 

The  foregoing  conclusions  are  based  on  the  assumption  that  each  subarray 
is  beamed  before  the  entire  array  is  beamed.  If  infinite-velocity  subarray 
beams  are  formed,  the  results  can  be  expected  to  differ  from  those  in  this 
report.  The  effects  of  infinite-velocity  subarray  beams,  and  number  of 
sensors  per  subarray,  are  currently  unde’-  investigation. 
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